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This  training  document  provides  guidance  on  the  development  of  frequency 
curves  from  annual  peak  discharges  that  are  segregated  into  two  populations. 
While  the  procedures  contained  in  this  document  use  annual  peaks  caused  by 
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CHAPTER  1.  INTRODUCTION 


The  development  of  frequency  curves  at  certain  locations  may  require 
special  treatment  when  the  events  are  caused  by  different  types  of  hydrologic 
phenomena  and/or  the  frequency  curve  exhibits  a  sudden  change  in  curvature. 

One  example  where  mixed  populations  occur  is  along  the  Atlantic  and  Gulf  coasts. 

In  these  areas  events  are  caused  by  intense  tropical  and  cyclonic  storms,  which 
are  referred  to  as  hurricane  and  non-hurricane  events.  Another  example  is  in 
the  Sierra  Nevada  region  of  California  where  rainfloods  tend  to  occur  November 
through  March  and  snowmelt  floods  generally  occur  April  through  July.  This 
training  document  discusses  the  development  of  a  frequency  curve  from  two  or 
more  sets  of  data  (populations)  that  originate  from  separate  causal  factors. 

The  terminology  used  in  this  document  is  as  follows.  When  the  frequency 
curve  is  derived  from  two  or  more  separate  frequency  curves,  each  developed 
from  a  separate  population,  the  resultant  curve  is  referred  to  as  the  combined- 
population  frequency  curve.  When  the  resultant  frequency  curve  is  derived 
directly  from  annual  peak  data  that  have  not  been  segregated  according  to 
causal  factors,  it  is  referred  to  as  a  mixed-population  frequency  curve. 

This  document  discusses  when  and  how  to  develop  a  combined-population  frequency 
curve  from  hurricane  and  non-hurricane  populations.  The  equations  can  be  used 
to  develop  e  combined-population  frequency  curve  from  other  mixed  populations. 

Chapter  2  discusses  the  merits  of  mixed-population  versus  combined-population 
frequency  analyses.  The  procedure  for  developing  a  combined-population  frequency 
curve  is  described  in  Chapter  3.  Examples  of  mixed-population  and  combined- 
population  frequency  curves  are  contained  in  Chapter  4. 
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CHAPTER  2.  WHEN  TO  USE  A  COMBINED-POPULATION  FREQUENCY  ANALYSIS 


The  combined-population  frequency  approach  should  be  considered  when  the 
frequency  curves  derived  from  mixed  populations  exhibit  rather  sudden  breaks  in 
the  curvature  of  the  frequency  curves.  Sometimes  unusually  large  or  small  skew 
coefficients  may  be  an  indication  of  mixed  populations.  Unusual  regional  skew 
coefficients  are  generally  considered  to  be  greater  than  0.7  and  less  than  -0.4. 

The  sudden  break  in  a  frequency  curve  is  often  caused  by  several  large  events 
that  depart  significantly  from  the  trend  of  the  rest  of  the  data.  These  large 
events  are  frequently  produced  by  a  different  type  of  hydrologic  phenomena; 
such  as  hurricanes  in  a  normally  rainfall  series,  rainflood  events  in  a  basically 
snowmelt  series  (U.S.  Geological  Survey,  1978),  or  thunderstorm  events  in  a 
basically  winter  rainstorm  series.  A  combined-population  analysis  is  often 
used  to  solve  this  problem,  but  because  of  the  additional  effort  required  to 
use  this  approach,  it  is  not  always  advantageous  to  do  so.  This  chapter 
discusses  the  considerations  involved  with  making  such  a  decision. 

The  primary  motivation  behind  a  combined-population  analysis  is  to  provide 
a  better  fit  between  the  analytically  derived  distribution  and  the  plotting 
positions  than  can  be  obtained  with  a  mixed-population  frequency  analysis.  If 
the  extreme  flood  events  are  considered  to  be  the  largest  in  a  time  period 
greater  than  the  systematic  annual  peak  flows,  then  procedures  contained  in 
Bulletin  17B,  Guidelines  for  Determining  Flood  Flow  Frequencies,  (U.S.  Water 
Resources  Council,  1981)  may  be  applicable.  The  weighting  of  high  events, 
based  on  an  extended  (histroic)perlod  of  record,  can  reduce  the  departure  of 
the  high  events  from  the  analytical  frequency  curve.  If  historical  information 
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Is  available  and  the  incorporation  of  these  data  in  a  frequency  analysis  provides 
a  good  fit  to  the  plotted  data,  then  a  combined-population  frequency  analysis 
may  not  be  warranted. 

When  the  historical  adjustment  does  not  provide  a  reasonable  fit,  or  if 
historical  information  is  not  available,  then  the  combined-population  frequency 
approach  should  be  considered.  If  it  is  not  clear  that  the  one  population  is 
responsible  for  the  sudden  change  in  curvature  in  a  fairly  large  number  of  cases, 
then  a  standard  frequency  analysis  using  the  mixed-population  approach  is 
preferred  for  three  reasons.  First,  it  may  be  difficult  to  identify  all  the 
events  for  each  population.  Second,  if  there  are  a  small  number  of  occurrences 
of  one  population,  the  resultant  frequency  curves  are  not  reliable  and  smoothing 
of  the  computed  statistics  is  required.  And  third,  much  effort  must  be  expended 
in  deriving  generalized  skew  coefficients  for  each  of  the  separate  populations. 

A  special  consideration  for  analysis  of  hurricane  and  non-hurricane  events 
is  the  size  of  the  drainage  area.  In  a  small  drainage  area  the  rainfall  intensity 
of  a  non-hurricane  event  can  often  be  equal  to  that  of  the  hurricane  event.  As 
the  size  of  the  basin  of  interest  increases,  the  chance  of  a  non-hurricane  event 
equaling  the  intensity  of  a  hurricane  event  decreases.  Therefore,  the  effect 
of  the  hurricane  events  on  the  small  drainage  areas  is  not  as  pronounced. 

The  Hydrologic  Engineering  Center  has  found  that  catchments  less  than  500  square 
miles  generally  will  not  require  a  special  hurricane  analysis.  The  drainage 
area  is  not  a  consideration  in  the  decision  to  segregate  rainfall  and  snowmelt 
events . 
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In  cases  where  the  sudden  departures  in  curvature  are  noted  in  some  stations 
but  not  in  others,  the  region  may  need  to  be  subdivided  into  two  separate  areas 
and  separate  regional  analyses  employed  in  each  area.  Care  must  be  taken  to  be 
sure  that  there  are  sufficient  stations  in  each  area  to  perform  a  regional 
analysis. 

Another  important  consideration  is  the  independence  of  events.  If  the  data 
in  one  of  the  series  is  not  independent  of  data  in  the  other  series,  then  a 
coincident  frequency  analysis  rather  than  a  combined-population  frequency  analysis 
will  be  warranted. 
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CHAPTER  3.  PROCEDURE  FOR  PERFORMING  A  COMBINED- 
POPULATION  FREQUENCY  ANALYSIS 


This  chapter  describes  the  procedure  for  developing  a  combined-population 
frequency  curve  from  hurricane  and  non-hurricane  Induced  events.  Section  3.1 
discusses  the  selection  of  events.  Section  3.2  summarizes  the  procedures  and 
provides  references  for  performing  a  standard  frequency  analysis.  Section 
3.3  describes  several  methods  for  determining  annual  frequency  curves  from  a 
set  of  events  that  do  not  occur  every  year.  The  procedure  for  combining  two 
frequency  curves  Is  reviewed  In  Section  3.4  and  the  development  of  regional 
relationships  used  to  develop  frequency  curves  at  ungaged  sites  is  mentioned 
in  Section  3.5.  Procedures  to  calculate  an  approximate  expected  probability 
adjustment  and  estimates  of  confidence  limits  are  contained  in  Sections  3.6 
and  3.7,  respecti vely. 


Section  3.1.  Data  Selection 

The  first  step  Is  to  obtain  the  necessary  data  to  perform  the  mixed- 
population  frequency  analysis.  Usually  the  annual  peaks  can  be  obtained 
directly  from  the  U.S.  Geological  Survey  (USGS)  Water  Supply  Papers.  Because 
of  the  effort  required  to  gather  data  for  the  combined-population  frequency 
analysis,  the  mixed-population  approach  should  be  completed  first  to  determine 
if  additional  analyses  are  warranted. 

The  collection  of  data  for  the  combined-population  frequency  analysis  Is 
the  next  step  and  Involves  the  determination  of  the  causal  factors  of  the 
events,  anff  the  Identification  of  the  largest  annual  event  In  each  population. 


Hurricane  events  can  be  Identified  by  studying  publications,  such  as  the 
U.S.  Department  of  Commerce's  report  on  Tropical  Cyclones  of  the  North  Atlantic 
Ocean  (1965),  that  specify  the  hurricane  tracks.  Intensity,  and  dates  of  occur¬ 
rence  at  selected  locations.  Next,  the  dates  of  these  events  can  be  compared 
with  flood  events  in  the  USGS  Water  Supply  Papers  to  determine  if  the  dates  of 
the  discharge  events  correspond  to  those  of  the  hurricanes.  This  methodology 
is  approximate  because  the  exact  location  and  the  areal  extent  of  the  hurricane 
is  not  known,  and  it  is  difficult  to  distinguish  between  a  remnant  of  a  hurricane 
event  and  a  non-hurricane  event.  Hurricane  events  do  not  occur  every  year  over 
most  drainage  areas  and  require  special  procedures  to  develop  a  frequency  curve 
as  discussed  in  Section  3.3. 

Another  typical  application  of  the  combined-population  frequency  analysis 
Is  the  division  of  the  year  into  seasons  or  months.  Bulletin  17B  states  that 
"separation  by  calendar  periods  In  lieu  of  separation  by  events  Is  not  consid¬ 
ered  hydrologlcally  reasonable  unless  the  events  In  the  separated  periods  are 
clearly  caused  by  different  hydrometeorological  conditions."  The  HEC  has  found 
that  if  the  data  are  segregated  Into  too  many  seasons,  then  one  or  more  of  the 
seasonal  frequency  curves  may  contain  one  or  two  large  events  and  many  small 
ones.  This  causes  the  seasonal  curves  to  have  a  very  steep  slope;  and  when  the 
seasonal  curves  are  combined  Into  a  single  annual  curve.  It  causes  the  upper 
end  of  the  annual  frequency  curve  to  be  unreasonably  high.  In  addition,  as 
stated  in  the  U.S.  Army  Corps  of  Engineers  Engineering  Manual  on  Hydrologic 
Frequency  Analysis  (1980).  "the  combined  curve  will  very  likely  fit  the  annual 
curve  only  in  the  middle  parts  of  the  curve,  and  the  lower  end  of  the  curve 
will  have  a  partial  duration  shape  as  many  small  events  have  been  included 


in  the  analysis." 


Section  3.2.  Standard  Frequency  Analysis 

If  a  complete  series  of  annual  peaks  can  be  identified,  a  standard  analy¬ 
sis  can  be  performed  to  determine  the  annual  peak  discharge  frequency  curve. 
These  procedures  have  been  extensively  documented  in  numerous  publications 
(U.S.  Army  Corps  of  Engineers,  1962  and  1975a;  and  U.S.  Water  Resources  Coun¬ 
cil,  1981),  and  are  briefly  summarized  in  this  section.  A  frequency  curve  is 
developed  from  the  annual  series  of  data  for  each  population  using  the  proce¬ 
dures  described  below. 


The  first  step  is  to  determine  graphical  plotting  positions  that  define 
the  exceedance  probability  associated  with  each  discharge.  The  annual  peak 
data  are  ranked  in  descending  order  and  a  plotting  position  is  determined 
using  one  of  several  different  equations.  One  of  the  most  common  is  the 
Weibull  plotting  position  equation  shown  below: 

p-m  ^ 

where  P  =  exceedance  probability  corresponding  to  the  event  of  rank  m 
m  =  rank  of  the  event 
N  =  number  of  events 


This  equation  was  developed  so  that  the  exceedance  probability  associated  with 
the  highest  ranked  event  would  be  correct,  on  the  average.  Equation  3.2  is 
another  commonly  used  plotting  position  which  is  an  approximation  of  the  Beard 
or  median  plotting  position  (U.S.  Arny  Corps  of  Engineers,  1962). 


n  m-0.3 
K  "  S+P 


(3.2) 
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The  median  plotting  position  was  developed  so  that  the  exceedance  probability 
associated  with  the  largest  event  would  have  an  equal  chance  of  being  too  high 
or  too  low.  Once  the  plotting  positions  have  been  determined,  the  exceedance 
probability  and  discharge  coordinates  are  plotted  on  the  appropriate  probability 
paper. 

An  analytical  frequency  curve  is  then  calculated  using  the  recommended 
probability  distribution.  The  U.S.  Water  Resources  Council  (1981)  recommends 
that  the  log-Pearson  type  III  distribution,  with  a  weighted  skew  coefficient, 
be  used  to  model  annual  peak  discharges.  However,  WRC's  conclusions  and 
generalized  skew  coefficient  map  were  based  on  annual  peak  data  that  were  not 
segregated  according  to  causal  factors.  If  the  log-Pearson  type  III  distri¬ 
bution  is  desired  to  model  a  segregated  series,  then  the  investigator  will  either 
need  to  accept  the  fundamental  uncertainty  of  a  calculated  skew  coefficient, 
or  perform  the  necessary  studies  for  developing  a  generalized  skew  relationship 
for  each  type  of  series.  Unless  the  annual  series  in  a  number  of  stations 
clearly  contain  non-zero  skew  coefficients,  a  log-normal  distribution  is  rec¬ 
ommended. 

The  analytical  frequency  curve  for  each  population  is  calculated  and 
plotted  along  with  the  corresponding  graphical  plotting  points.  The  expected 
probability  adjustment  and  the  confidence  limits  for  the  analytical  curve  are 
not  determined  until  the  combined-population  frequency  curve  has  been  derived, 
as  described  in  Sections  3.6  and  3.7. 


Section  3.3.  Development  of  Frequency  Curves  From  a  Truncated  Series 


Special  frequency  analyses  are  required  when  events  In  a  series  do  not 
occur  every  year.  This  section  discusses  two  procedures  that  have  been  used  to 
develop  hurricane  frequency  carves  In  several  HEC  studies. 

The  first  procedure  used  at  the  HEC  was  developed  by  Beard  (U.S.  Army  Corps 
of  Engineers,  1958).  He  recommended  that  a  standard  frequency  analysis  be  per¬ 
formed  on  the  hurricane  events  yielding  a  curve  based  on  the  number  of  hurricane 
events  per  100  events.  This  curve  Is  a  conditional  frequency  curve.  Identical 
in  concept  to  the  one  discussed  by  the  U.S.  Water  Resources  Council  (1981). 

The  exceedance  probabilities  of  this  curve  are  then  multiplied  by  N^/Nj,  where 
is  the  number  of  hurricane  events  and  Ny  is  the  total  number  of  years  of 
record.  While  this  adjustment  can  dramatically  affect  the  lower  end  of  the 
hurricane  curve,  it  causes  only  a  moderate  shift  of  the  frequency  curve  at  its 
upper  end.  This  technique  has  not  been  used  in  many  applications  because  it  is 
usually  considered  valid  when  only  less  than  25  percent  of  the  data  is  missing. 
Because  this  is  seldom  the  case  with  hurricane  events,  an  alternative  procedure 
has  been  used  by  the  HEC  in  several  applications  (U.S.  Army  Corps  of  Engineers, 
1965  and  1975b). 

The  first  step  In  this  alternative  procedure  Is  to  compute  plotting  posi¬ 
tions  of  the  data  series  using  either  Equation  3.1  or  3.2  In  the  same  manner  as 
described  In  Section  3.2  except  that  N  Is  the  number  of  years  rather  than  the 
number  of  events.  The  frequency  curve  Is  then  developed  by  drawing  a  best-fit 
line  through  the  plotting  positions.  This  line  can  be  based  either  by  eye  or  by 
a  modified  regression  technique  which  provides  a  more  rigorous  mathematical 
estimate  (U.S.  Army  Corps  of  Engineers,  1959).  However,  the  modified  regression 
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procedure  places  equal  weight  on  each  of  the  data  points,  so  that  one  outlier  can 
drastically  affect  the  derived  line.  The  slope  of  the  line  developed  by  one  of 
these  procedures  is  the  standard  deviation.  The  mean  of  the  hurricane  events 
is  obtained  by  extending  the  adopted  line  and  noting  the  discharge  associated 
with  the  0.50  exceedance  probability. 

Due  to  the  small  sample  typically  used  to  develop  a  frequency  curve  in 
this  manner,  there  is  a  great  deal  of  uncertainty  in  both  the  mean  and  standard 
deviation.  This  deviation  may  vary  considerably  at  different  geographic  locations. 
Therefore  the  mean  and  standard  deviation  are  often  plotted  versus  the  drainage 
area  at  each  gaged  site  to  provide  a  basis  for  selecting  a  regional  value. 

Different  mean  and  standard  deviation  relationships  may  be  adopted  for  different 
river  systems  (U.S.  Army  Corps  of  Engineers,  1975b)  or  a  single  relationship  for 
the  mean  and  a  single  value  for  the  standard  deviation  may  be  adopted  for  the 
entire  region  (U.S.  Army  Corps  of  Engineers,  1965).  A  zero  skew  coefficient 
is  generally  adopted  unless  there  is  a  regional  trend  that  can  be  rationalized 
to  be  caused  by  known  climatic  or  basin  characteristics. 

Because  of  the  difficulty  in  identifying  hurricane  events,  and  because  of 
their  small  sample  size,  there  is  a  great  deal  of  uncertainty  in  estimating  the 
mean  and  standard  deviation  at  individual  sites.  Therefore,  even  though  the 
smoothing  techniques  may  be  highly  subjective,  they  are  desirable. 

Section  3.4.  Combining  Frequency  Curves 

The  procedure  for  combining  frequency  curves  developed  from  independent 
annual  series  has  been  widely  documented  (U.S.  Army  Corps  of  Engineers,  1958, 

1965,  1975b,  and  1980).  The  general  equation  for  combining  multiple  frequency 
curves  is: 

\  ■  l-(l-P,)(l-P2)...(l-P11)  ■  l-MI-P,)  ,3.3, 
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where  Pc  is  the  exceedance  probability  of  the  combined-population  frequency 
curve  for  the  selected  discharge 

P-j,P2,...Pn  are  the  exceedance  probabilities  associated  with  a  selected 
discharge  from  frequency  curve  numbers  1,2,  through  n 
n  is  the  number  of  frequency  curves  that  are  combined 


If  only  two  curves  are  combined,  then  Equation  3.3  reduces  to: 


Pc  -  VW, 


(3.4) 


Equations  3.3  and  3.4  are  only  valid  when  each  of  the  frequency  curves  used  to 
develop  a  combined  curve  are  assumed  to  be  independent.  (See  Appendix  I  for 
development  of  these  equations.) 


Section  3.5.  Combined-Population  Frequency  Curves  in  Ungaged  Areas 

This  section  describes  several  procedures  for  developing  frequency  curves 
at  ungaged  sites  using  combined-population  frequency  analysis  results  at  gaged 
sites.  The  investigator  is  faced  with  the  choice  of  determining  regional 
relationships  that  calculate  a  mixed-population  curve  directly  or  developing 
regional  relationships  that  calculate  a  frequency  curve  for  each  population. 
While  the  latter  procedure  is  more  theoretically  appealing,  there  is  often 
greater  uncertainty  in  individual  frequency  curves  developed  from  the  separate 
populations,  and  consequently  the  regional  relationships  are  apt  to  contain  a 
great  deal  of  uncertainty.  Unless  the  investigator  is  very  confident  in  the 
analyses  and  data  used  to  develop  each  of  the  separate-population  frequency 
curves,  a  mixed-population  approach  might  be  warranted  because  it  is  simpler 
to  perform  and  may  provide  the  same  level  of  reliability. 
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The  HEC  developed  separate  regional  relationships  for  hurricane  and  non¬ 
hurricane  events  on  the  Kanawha  River  (U.S.  Army  Corps  of  Engineers,  1965). 

Three  empirical  curves  were  developed  that  described  the  drainage  area  versus 
the  mean  of  the  logarithm  of  the  hurricane  events,  the  mean  of  the  logarithm  of 
the  non-hurricane  events,  and  the  standard  deviation  of  the  logarithm  of  the 
non-hurricane  events.  Because  there  were  inadequate  data  on  hurricane  floods, 
a  uniform  standard  deviation  of  the  logarithms  of  the  hurricane  events  was 
adopted.  These  relationships  could  be  developed  for  hurricane  and  non-hurricane 
frequency  curves  and  then  be  combined  using  Equation  3.4. 

An  alternative  approach  was  used  in  the  Tropical  Storm  Agnes  Study  (U.S. 

Army  Corps  of  Engineers,  1975b).  In  this  report,  separate  regional  relationships 
for  calculating  hurricane  and  non-hurricane  events  were  developed  for  ungaged 
areas  along  the  major  rivers.  Regression  equations  were  determined  from  the 
mean  and  standard  deviation  of  the  logarithm  of  the  mixed-population  frequency 
curves.  These  equations,  along  with  a  regional  skew  map,  were  developed  for  use 
in  the  ungaged  areas  that  were  not  along  the  major  rivers.  By  developing  a  mixed- 
population  curve  directly  at  the  ungaged  sites,  the  tremendous  uncertainty 
involved  with  regionalization  of  the  hurricane  events  was  avoided.  However, 
the  use  of  this  approach  did  not  address  the  problem  of  the  sharply  skewed 
frequency  curve. 


Section  3.6.  Expected  Probability 

The  expected-probability  methodology,  as  proposed  by  Beard  (U.S.  Army  Corps 
of  Engineers,  1962)  adjusts  frequency  estimates  so  that  the  average  of  the 
exceedance  probabilities  for  many  different  sites  is  closer  to  the  true  population 
exceedance  probability.  It  was  developed  assuming  a  normal  distribution  but  has 
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also  been  used  to  adjust  frequency  curves  developed  using  a  log-Pearson  type  III 
distribution.  Hardison  and  Jennings  (1972)  have  shown  that  this  adjustment 
reduces  the  bias  In  samples  drawn  from  a  log-Pearson  type  III  distribution  with 
a  known  skew  coefficient.  Lloyd  (1978)  indicated  that  similar  estimating  proce¬ 
dures  could  be  developed  for  other  distributions.  However,  he  also  indicated 
that  their  mathematical  derivation  is  extremely  complex. 

The  underlying  distribution  of  the  combined-population  frequency  curve 
developed  using  Equation  3.4  is  unknown.  Therefore  the  correct  adjustment  to 
cause  the  average  of  the  probabilities  from  many  sample  frequency  curves  to 
equal  the  "true"  population  exceedance  probability  is  also  unknown.  The  trade¬ 
off  becomes  one  of  whether  a  possibly  incorrect  mathematical  adjustment  should 
be  applied,  or  no  adjustment  at  all.  The  HEC  has  made  the  adjustment  in  their 
studies  because  even  though  the  adjustment  may  not  be  known  with  much  uncertainty, 
it  was  felt  that  it  would  cause  less  bias  in  frequency  estimates  than  no  adjust¬ 
ment. 

The  expected  probability  adjustment  can  be  made  using  Chart  40  in  "Statis¬ 
tical  Methods  in  Hydrology"  (U.S.  Artny  Corps  of  Engineers,  1962)  or  by  employing 
Equation  11-1  in  the  Water  Resource  Council  Guidelines  (1981)  with  a  cumulative 
distribution  of  the  student's  t  distribution  (Benjamin  and  Cornell,  1970).  Both 
procedures  for  making  the  expected  probability  adjustment  require  using  the  num¬ 
ber  of  events,  n,  in  the  sample.  For  a  combined-population  frequency  curve  at 
a  gaged  site,  the  HEC  has  defined  n  as  the  larger  of  the  number  of  events  used 
to  develop  the  non-hurricane  or  the  hurricane  frequency  curve.  However,  at  un¬ 
gaged  sites,  the  mean  and  standard  deviation  are  often  developed  from  separate 
regional  equations.  The  worth  of  the  mean  and  standard  deviation  will  depend 
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upon  the  correlation  coefficient  and  the  number  of  years  of  record  at  the  gaged 
sites  used  to  develop  these  relationships.  The  HEC  has  been  unable  to  develop 
a  theoretically  appealing  mathematical  relationship  to  establish  an  equivalent 
period  of  record  at  an  ungaged  site.  Typically,  the  average  number  of  years  of 
record  at  the  gaged  sites  in  the  region  has  been  calculated  and  used  to  make  the 
expected  probability  adjustment  at  ungaged  locations. 

Section  3.7.  Confidence  Limits  for  Combined-Population  Frequency  Curves 

Confidence  limits  for  the  normal  distribution  can  be  calculated  using  a 
non-central  t  distribution  (Resnikoff  and  Lieberman,  1957).  The  non-central  t 
has  also  been  employed  to  calculate  confidence  limits  for  a  log-Pearson  type  III 
distribution  with  a  skew  coefficient  between  ±0.5  (U.S.  Water  Resources  Council, 
1981).  However,  in  the  latter  case,  only  the  uncertainties  in  the  mean  and  stan¬ 
dard  deviation  are  accounted  for.  When  the  underlying  distribution  is  not  normal, 
the  calculation  of  confidence  limits  based  on  the  non-central  t  may  not  be  theo¬ 
retically  valid. 

Establishing  confidence  limits  for  a  combined-population  frequency  curve 
has  the  same  theoretical  difficulties  as  making  the  expected  probability  adjust¬ 
ment:  either  calculate  a  possibly  incorrect  mathematical  confidence  limit  or 
none  at  all.  The  HEC  has  used  two  alternate  procedures  to  calculate  approximate 
confidence  limits  for  the  combined-population  frequency  curve  as  described  in 
this  document.  These  techniques  are  only  suggested  solutions  to  the  problem  at 
this  time.  Further  analysis  needs  to  be  carried  out  to  verify  these  procedures 
or  develop  new,  more  reliable  procedures. 
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Method  1 


The  first  procedure  for  calculating  confidence  limits  employs  Exhibit  6 
in  "Statistical  Methods  in  Hydrology"  (U.S.  Army  Corps  of  Engineers,  1962). 
This  chart  tabulates  an  error  of  the  estimated  value,  EM  D  which  is  a 
function  of  the  years  of  record  N,  the  exceedance  probability  P,  and  the 
confidence  level  c.  The  confidence  limits  can  then  be  calculated  using  the 
equation: 

XP,C  '  XP  +  EN.P.C  *  S  <3-5> 

where  S  is  the  standard  deviation  of  the  frequency  curve,  Xp  is  the  logarithm 
of  the  discharge  at  the  exceedance  probability  P,  and  XD  is  the  logarithm  of 
the  confidence  limit  with  exceedance  probability  P.  The  tabulated  values  of 
EN  P  c  are  P°sitive  for  the  uPPer  confidence  limits  and  negative  for  the  lower 
limits.  If  values  of  EN  p  c  are  needed  which  are  not  tabulated  in  Exhibit  6 
(U.S.  Army  Corps  of  Engineers,  1962),  they  can  be  calculated  by  the  following 
equations: 


v  *  Vp.c 

P,C  /ITT  (3.7) 

where  Kp  is  the  normal  deviate  for  exceedance  probability  P,  and  tN  p  c  is  the 
non-central  t  value.  Values  of  xp  a  non-central  t  argument,  for  various  degrees 
of  freedom  (N-l  in  this  application)  can  be  found  in  Resnikoff  and  lieberman  (1957), 
and  Kp  is  found  in  numerous  sources  (e.g.,  U.S.  Water  Resources  Council,  1981). 

If  non-central  t  tables  are  not  available.  Equations  9-4  through  9-6  in  the 
Water  Resources  Guidelines  can  be  used  to  determine  approximate  values  of  Kp  . 


where  Sc  Is  the  standard  deviation  of  the  combined-population  frequency  curve 
for  the  discharge  Q;  Ph^q  is  the  exceedance  probability  from  the  hurricane 
frequency  curve  associated  with  the  discharge  Q;  P^  q  Is  the  exceedance 
probability  from  the  non-hurricane  frequency  curve  associated  with  discharge 
Q;  SH  Is  the  standard  deviation  of  the  hurricane  frequency  curve,  and  Is 
the  standard  deviation  of  the  non-hurricane  frequency  curve.  As  shown  In 
the  example  In  Chapter  4,  this  procedure  will  not  always  yield  valid  confidence 
limits.  If  there  Is  a  large  difference  between  the  hurricane  and  non-hurricane 
standard  deviations  and  skew  coefficients,  this  can  lead  to  an  irregularly 
shaped  lower  confidence  limit  below  the  intersection  of  the  two  curves. 

Method  2 


An  alternate  procedure  1$  to  calculate  separate  confidence  limits  for  the 
hurricane  and  non-hurricane  frequency  curves.  These  limits  can  be  combined 
using  equation  3.4.  While  this  method  avoids  the  problems  of  the  Irregularly 
shaped  confidence  limits.  It  does  yield  confidence  bounds  that  are  perhaps  too 
close  to  the  combined-population  frequency  curve. 

If  the  standard  deviation  and  skew  coefficients  are  reasonably  close,  the 
method  using  Equations  3.5  -3.12  should  provide  satisfactory  results.  When 
this  Is  not  true  the  second  method  that  uses  Equation  3.4  to  combine  the 
separately  derived  confidence  limits  may  provide  more  consistent  results.  In 
either  case  the  investigator  must  carefully  examine  the  derived  confidence 
limits  to  be  sure  they  appear  reasonable  and  remember  that  both  procedures 
provide  only  rough  estimates  of  the  actual  uncertainty  in  the  combined- 
population  frequency  curve. 
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CHAPTER  4.  EXAMPLES  OF  MIXED-POPULATION  AND  COMBINED-POPULATION  FREQUENCY 
ANALYSIS 

This  chapter  contains  examples  of  mixed-population  and  combined-population 
frequency  analysis  for  a  stream  gaping  station  on  West  Conewago  Creek  near 
Manchester,  Pennsylvania.  Section  4.1  summarizes  the  data  and  Sections  4.2 
and  4.3  describe  the  computations  for  mixed  and  combined  analyses,  respectively. 
The  frequency  curves  derived  In  these  Chapters  are  for  Illustrative  purposes 
only;  they  are  not  to  be  construed  as  the  recommended  curves  Tor  West  Conewago 
Creek. 


Section  4.1.  Data  Selection 

Figure  4.1  contains  an  excerpt  from  a  U.S.G.S.  water  supply  paper  and 
a  list  of  additional  annual  peaks  for  West  Conewago  Creek  near  Manchester,  PA 
for  the  years  1929  to  1972.  A  separate  analysis  has  Indicated  that  five 
hurricane  events  were  interspersed  throughout  the  44  years  of  recorded  data. 
Both  the  hurricane  and  non-hurrl cane  annual  peak  discharges  for  these  years 
are  shown  In  Figure  4.2. 


18 


5740 .  Vest  Cor.ewago  Creek  rear  Manchester,  Pa, 

(Published  as  "Conewago  Creek'1  prior  to  1932) 

Location .--Lat  40"04'551',  long  76 “43 '10”,  500  ft  upstream  from  bridge  or.  State 
highway  24,  0.7  mile  downstream  from  Little  Conewagc  Creek, and  1.5  miles 
north  of  Manchester,  York  County. 

Drainage  area.— 510  eq  ml. 

Cage. "Recording.  Datum  of  gage  la  263.68  ft  above  mean  sea  level,  datum  of 
*1989. 

Stage-discharge  re lat Ion. --Defined  by  current-meter  measurements. 

Bar.kfull  stage.--?  ft. 

Remarks  .--Base  f'or  partlal-duratlor.  series,  10,600  cl's. 


Peak  stages  and  discharges  of  West  Cor.ewago  Creek  near  Manchester,  Pa. 


Water 

year 

Date 

Cage 

height 

(feet) 

Discharge 

(cfa) 

Water 

year 

— 

Date 

Gage 

height 

(feet) 

Discharge 

(cfe) 

.929 

Mar. 

6, 

1929 

11.70 

n,eoo 

i 

1945 

.July  19, 

1545 

11.39 

*  11,100 

Apr. 

17, 

1929 

15.31 

20,300 

May 

3, 

1929 

13.79 

16,500 

1! 

1946 

Nov.  29, 

1945 

15.86 

21,600 

June  2, 

1546 

15.74 

21,000 

193C 

Oct. 

2, 

1929 

12.56 

13,  '700‘ 

"s 

Mar. 

8, 

1930 

11.16 

10,600 

1947 

May  22 , 

2947 

13.66 

16,000 

1931 

Apr. 

2, 

1931 

9.14 

6,850 

li 

ji 

1945 

Jan.  2 , 

1946 

10.62 

9,980 

1932 

Mar. 

28, 

1932 

12.12 

11,900 

i 

1949 

Dec.  30, 

:®4e 

13.74 

16,000 

Jar. .  6 , 

1949 

12.96 

14,400 

1533 

Oct. 

13, 

1932 

13.37 

14,200 

Nov. 

1, 

1932 

12.62 

12,500 

1550 

Mar.  22, 

1950 

22.51 

23,300 

'  Apr. 

20, 

1933 

13.26 

14, IOC 

'1 

Kay  19, 

1950 

12.16 

12,700 

Aug. 

24, 

1833 

24.14 

47,600* 

!« 

1951 

Nov.  26, 

1950 

12.53 

13,300 

1534 

Sept. 

15, 

1934 

13.71 

15,300 

.■j 

Dec.  4, 

1950 

13.91 

16,400 

Sept. 

17. 

1934 

17.41 

24,900 

'j 

Feb.  7, 

1551 

11.67 

11,700 

Sept 

30, 

1934 

17.20 

24,400 

ii 

Feb.  21, 

1951 

11.64 

11,500 

1935 

Dec. 

2, 

1934 

I5.ee 

20,  700 

ij 

1952  * 

Pel*.  4, 

1552 

11.93 

12,100 

Kar.  11, 

1552 

13.88 

16,400 

1936 

Mar. 

12, 

1936 

13.02 

13,700 

ii 

Apr.  28, 

1952 

11.72 

11,700 

Mar. 

19, 

1936 

al7.oe 

- 

Apr. 

6, 

1936 

12.53 

13,500 

j' 

.553 

Nov.  22, 

1552 

13.96 

16 , 7QC 

June 

13, 

1936 

11.93 

11.40C 

ii 

Jan.  24, 

1553 

12.66 

13,700 

1537 

Feb. 

22, 

1937 

11.73 

12,100 

!i 

1954 

Mar.  2, 

1954 

8.30 

5,740 

Apr. 

27, 

1937 

12.06 

12,900 

/! 

1955 

Mar.  22, 

1955 

14.10 

It.SOC 

193c 

Oct. 

23, 

1937 

11,27 

n,20o 

ii 

Aug.  13, 

1955 

11.29 

10,900 

Nov. 

13, 

1937 

13.82 

16,800 

1556 

Oct.  14, 

1955 

12.92 

24,200 

1939 

Feb. 

4, 

1939 

13.70 

16,500 

Mar. 

1, 

1939 

11.16 

11,000 

1957 

Dec.  15, 

1956 

11.37 

11,100 

! 

Apr .  6 , 

1957 

11.34 

10,900 

1940 

Apr. 

9, 

1940 

12.16 

12,500 

Apr. 

20, 

1940 

15.95 

21,300 

) 

2956 

Dec.  21, 

2957 

23.77 

26, 20C 

Sept- 

1, 

1940 

11.63 

11,200*! 

Dec.  27, 

1557 

13.14 

14,600 

Feb.  2c, 

1958 

IF  97 

14,400 

1541 

Apr. 

e, 

1941 

11.16 

10,400 

ii 

Mar.  26, 

1956 

12.03 

12 , 300 

J  6. 

1953 

U.74 

11,700 

1942 

May 

22, 

1942 

12.79 

13,800 

H 

June 

S, 

1942 

12.59 

13,400 

n 

195. 

Jan.  22, 

1959 

10.07 

8,720 

Aug. 

18. 

1542 

14.26 

1 “ , 400 

I960 

Apr.  4, 

1560 

12.55 

lS.500 

1543 

Dec . 

30, 

1942 

14.05 

ie,9rc 

'! 

1961 

Feb.  26, 

1561 

13.00 

14,400 

1944 

Nov. 

9, 

1943 

17.33 

2  5,500 

Apr.  13, 

1961 

12.84 

14,000 

Mar. 

13, 

1944 

11.36 

10 , 60C 

' 

Mar. 

24, 

1944 

11.39 

10,900 

May 

7, 

1944 

Lil-.l7  j 

10 , 800 

a  Backwater  from  Ice. 


1962 

11000 

1968 

16400 

1963 

15000 

1969 

15500 

1964 

11500 

1970 

21300 

1965 

18100 

1971 

15700 

1966 

16000 

1972 

81700* 

1967 

19000 

*  Hurricane  events 

Figure  4.1  Annual  peak  discharges  at  West  Conewago  Creek 

(U.S.  Geological  Survey,  1968) 
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Year 

Non-hurricane  Discharge 
(cfs) 

Hurricane  Discharge 
(cfs) 

1930 

10800 

13700 

1933 

14100 

47600 

1940 

21300 

11200 

1955 

16900 

10900 

197? 

12700 

81700 

Figure  4.2 

Hurricane  and  non-hurricane 

annual  peak  discharges 

Section  4.2.  Mixed-Population  Frequency  Analysis 

The  annual  peaks  for  West  Conewago  Creek,  irrespective  of  causal  factors, 
were  provided  as  input  to  the  Flood  Flow  Frequency  Computer  Program  (U.S.  Army 
Corps  of  Engineers,  1976)  which  was  used  to  perform  a  frequency  analysis.  A 
generalized  skew  coefficient  of  0.5  was  obtained  from  the  map  in  Bulletin  17B 
(U.S.  Water  Resources  Council,  1981).  The  plotting  positions  are  shown  in 
Figure  4.3,  the  final  statistics  and  frequency  ordinates  in  Figure  4 .  - .  anu 
frequency  curve  in  Figure  4.5.  It  is  evident  that  the  two  largest  discharges, 
which  are  both  hurricane  events,  caused  a  very  high  calculated  skew  coefficient. 
Both  these  events  depart  significantly  from  the  analytical  curve.  Either  his¬ 
torical  information  should  be  sought  to  determine  if,  in  fact,  the  hurricane 
events  were  the  largest  over  a  historical  period  greater  than  the  44  years  of 
recorded  data,  or  a  combined-population  analysis  should  be  performed. 

Because  there  was  no  readily  accessible  historical  information  at  West 
Conewago  Creek,  information  at  adjacent  gages  was  sought.  The  closest  long 
term  stream  was  located  on  the  Susquehanna  River  hear  Harrisburg,  Pennsylvania, 
and  had  a  continuous  record  since  1889. 
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Figure  4.3  Nixed  frequency  analysis  for  systematic  data  at  West  Conewago  Creek 
-  plotting  positions 

Illustrative  Example 
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Figure  4.4  Mixed  frequency  analysis  for  systematic  data  of  West  Conewago  Creek 
-  statistics  and  discharges  for  selected  exceedance  probabilities. 


FINAL  RESULTS 

-FREOUtAtr  FLU T  -  UEST  CONEUAGO  CREEK  HR  MANCHESTER, PA 
BASED  ON  CONFUTED  VALUES,  FLOW  IN  CUBIC  FEET  Fek  SECOND 


The  1972  Agnes  event  was  the  largest  in  the  period  of  record.  It  was  as¬ 
sumed  that  the  1972  Agnes  event  at  West  Conewago  Creek  was  also  the  largest 
since  1889.  The  Flood  Flow  Frequency  computer  program  was  run  using  this  his¬ 
torical  information  and  the  resultant  statistics  and  frequency  curve  are  shown 
in  Figures  4.6  ana  4.7.  In  this  particular  case,  tne  historical  adjustment 
lowers  the  upper  end  of  the  frequency  roughly  ten  percent. 

The  plotting  positions  of  the  two  highest  events  are  still  well  above  the 
analytical  frequency  curve.  The  following  section  illustrates  the  use  of  a 
combined-population  frequency  analysis  to  handle  this  problem. 


Section  4.3.  Combined-Population  Frequency  Analysis 

This  section  illustrates  the  basic  steps  used  to  develop  a  combined 
population  frequency  curve  using  the  data  shown  in  Figures  4.1  and  4.2. 


Section  4.3.1. 


nent  of  a  Non-Hurricane  Frequency 


The  non-hurricane  frequency  curve  was  derived  from  the  44  years  of  non¬ 
hurricane  events  using  the  Flood  Flow  Frequency  computer  Droqram.  The  plotting 
positions,  statistics  and  discharges  for  selected  exceedance  probabilities, 
and  the  frequency  curve  are  shown  in  Figures  4.8,  4.9,  and  4.10,  respectively. 


Section  4.3.2. 


Only  five  hurricane  events  were  noted  during  the  44  years  of  recorded 
flows,  therefore  the  hurricane  frequency  curve  will  be  based  on  a  truncated 
series.  The  frequency  curve  is  developed  using  the  second  procedure  discussed 
in  Section  3.3.  Plotting  point  positions,  tabulated  in  Figure  4.11,  were 
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Figure  4. 
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Mixed  frequency  analysis  for  systematic  and  historical  Information 
at  West  Conewago  Creek  -  statistics  and  discharges  for  selected 
exceedance  probabilities. 

Illustrative  Example 
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FI  ML  RESULTS 

-FREQUENCY  PLOT  -  WEST  CONEUAGO  CREEK  NR  MANCHESTER, PA 
RASED  ON  COMPUTED  VALUES,  FLOU  IN  CUBIC  FEET  PER  SECOND 


Figure  4.7  Mixed  frequency  analysis  for  systematic  and  historical  Information  at  West  Conewago  Creek 
-  plot  of  frequency  curve. 

Illustrative  Example 
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Figure  4.8  Non-hurricane  frequency  analysis  at  West  Conewago  Creek 
-  plotting  positions. 

Illustrative  Example 
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Figure  4.9  Non-hurricane  frequency  analysis  at  West  Conewago  Creek 

•  statistics  and  discharges  for  selected  exceedance  probabilities 


28 


Peak  flow  in  thousand 


Year 

Event 

Plotting  Position 

1972 

81700 

0.0222 

1933 

47600 

0.0444 

1930 

13700 

0.0667 

1940 

11200 

0.0889 

1955 

10900 

0.1111 

Figure  4.11  Hurricane  plotting  positions 


calculated  using  Equation  3.1.  N  was  defined  as  44,  the  number  of  years  of 
record.  These  plotting  points  are  drawn  on  Figure  4.12  and  an  eye-fit  curve  was 
drawn.  In  addition,  a  modified  regression  line  was  calculated  using  a  procedure 
recommended  by  Beard  (U.S.  Army  Corps  of  Engineers,  1959).  The  calculations  are 
shown  in  Figure  4.13  and  the  calculated  regression  line  is  also  drawn  on  Figure 
4.12.  The  eye-fit  curve  was  selected  for  inclusion  in  the  combined-frequency 
curve  analysis.  The  modified  regression  line  was  judged  to  be  too  steep  because 
of  the  following.  First,  it  is  quite  possible  that  the  hurricane  event  is  the 
largest  in  a  period  longer  than  the  44  recorded  years  of  data  as  discussed  in 
Section  4.2.  Second,  the  modified  regression  line  is  quite  sensitive  to  the 
magnitudes  of  the  smaller  hurricane  events.  Third,  other  regional  studies  have 
shown  that  for  basins  of  similar  size,  the  standard  deviation  varies  from  0.6  to 
1.2,  and  the  slope  of  the  eye-fit  line  falls  comfortably  in  the  middle  of  this 
range.  The  adopted  hurricane  frequency  curve  with  a  mean  of  2.9731  and  a  stan¬ 
dard  deviation  of  0.871  is  displayed  with  the  non-hurricane  curve  in  Figure  4.10. 

Section  4.3.3  Development  of  a  Combined-Population  Frequency _Curve 

The  hurricane  and  non-hurricane  frequency  curves  that  were  developed  in  the 
previous  two  sections  are  combined  using  Equation  3.4.  Discharges  are  selected 
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Peak  flow  In  thousand 


I.  Tabulation  of  log  Q  versus  the  corresponding  K  deviate. 


II. 


Event 

Q  cfs 

log  Q 
(X) 

K-normal  deviate  corresponding 
to  plotting  position 

81700 

4.9122 

2.0124 

47600 

4.6776 

1.7041 

13700 

4.1367 

1.5242 

11200 

4.0492 

1.3629 

10900 

4.0374 

1.2327 

For  this  example  X  equals  log  Q,  N  is  the  number  of  events,  and  K  can  be 
obtained  from  tables  of  the  standardized  normal  distribution  (Benjamin 
and  Cornell,  1970). 

Calculation  of  the 

standard  deviation  of  the  frequency  curve. 

/  EX2  - 

(rx)2/n 

SX  0.40506  _  n 

\l  Z Re¬ 

(rK)*/n 

SK  ~  0.30513  ‘ 

calculation  of  the 

mean  of  the  frequency  curve. 

Mx  =  IX/N 

=  4.3626 

Mk  =  ZK/N 

=  1.5673 

Mean  of  the  frequency  curve  = 

Mx  -  S  Mk  «  4.3626  -  1.328  (1.5673) 

— 

2.813 

IV.  Equation  of  frequency  curve  (regression  line). 
X  =  1.328  K  +  2.813 


Figure  4.13  Calculation  of  modified  regression  line 
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that  define  the  range  of  events  from  both  curves  and  the  corresponding 
exceedance  probabilities  are  picked  off  the  frequency  curve  and  used  to 
compute  the  exceedance  probability  of  the  combined  curve.  In  this  example  an 

additional  step  is  used  to  obtain  the  qreatest  possible  accuracy.  The  normal 
(Kp)  and  Pearson  Type  III  (Kg  p)  deviates  are  computed  using  Euqations  4.1 
and  4.2  for  the  hurricane  and  non-hurricane  curves,  respectively. 


(4.2) 

G  ,P  T 

where  and  ^  are  mean  of  the  logarithms  for  the  hurricane  and  non-hurricane 
events  respectively,  and  SH  and  are  the  standard  deviations  of  the  logarithms 
for  the  hurricane  and  non-hurricane  events  respectively.  The  exceedance 
probabilities  corresponding  to  these  deviates  can  be  found  using  tables  of  the 
normal  or  Pearson  type  III  distributions  (Benjamin  and  Cornell,  1970;  U.S.  Water 
Resources  Council,  1981"). 

Figure  4.14  summarizes  the  calculations  to  develop  a  combined-population 
frequency  curve  for  West  Conewago  Creek.  Column  1  contains  the  selected  discharges 
and  columns  2  and  4  contain  the  normal  and  Pearson  type  III  deviates  calculated 
using  Equations  4.1  and  4.2,  respectively.  The  corresponding  exceedance  prob¬ 
abilities  are  used  In  Equation  3.4  to  calculate  the  combined-population  exceed¬ 
ance  probability  In  column  6  that  corresponds  to  the  discharge  tabulated  In 
column  1. 
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Section  4.3.4.  Expected  Probability  Calculations 


The  expected  probability  adjustment  Is  calculated  for  the  combined-popula¬ 
tion  frequency  curve  using  44  years  of  record.  An  expanded  version  of  Exhibit  40 

in  "Statistical  Methods  in  Hydrology"  (U.S.  Army  Corps  of  Engineers,  1962)  was 

used  to  determine  the  P^  values  as  shown  below  In  Figure  4.15,  and  plotted  on 
Figure  4.16. 


p 

09 

PN 

Q 

0.0100 

0.0132 

99900 

0.0400 

0.0453 

57800 

0.0500 

0.0556 

26600 

0.1000 

0.1059 

22700 

0.3000 

0.3034 

18300 

0.4000 

0.4017 

16900 

0.5000 

0.5000 

15700 

0.6000 

0.5983 

14500 

0.7000 

0.6966 

13300 

0.9000 

0.8941 

9950 

0.9500 

0.9444 

8590 

Figure  4.15  Tabulation  of  P^  vs  ?m  for  the  combined-population  curve. 
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Section  4.3.5.  Confidence  Limits  Calculations 


The  two  methods  for  calculating  confidence  limits  that  were  described  in 
Section  3.7  are  used  in  this  section  to  develop  5  and  95  percent  confidence 
limits  for  West  Conewago  Creek. 


Method  1 

Equations  3.5,  3.8,  3.9  and  3.12  are  used  to  determine  the  confidence  limits 
as  illustrated  in  columns  7  through  21  in  Figure  4.14.  The  first  step  is  to  de¬ 
termine  equivalent  values  of  the  Pearson  type  III  deviates,  K^.  p  which  are 
needed  in  Equations  9-4  through  9-6  in  Bulletin  17B  to  calculate  values  of  Kp  ^ 
in  Equations  3.8  and  3.9.  z^  is  the  normal  deviate  that  is  associated  with  the 
level  of  significance  c. 

Values  of  Kp,  in  column  7  are  the  equivalent  normal  deviates  corresponding 
to  the  probabilities,  Pc,  in  column  6  of  the  combined  curve.  In  order  to  get 
the  equivalent  Pearson  deviate  of  the  combined  curve  KG  p  the  equivalent  skew 
coefficient  (Gc)  must  be  determined  for  each  discharge.  This  is  accomplished 
using  an  equation  identical  to  Equation  3.12,  except  the  skew  coefficient  of 
the  hurricane  and  non-hurricane  frequency  curves  are  used  in  lieu  of  the  stan¬ 
dard  deviations.  Column  8  is  the  relative  weight  of  the  non-hurricane  probabil¬ 
ities,  and  column  9  is  the  relative  weight  of  the  hurricane  probabilities. 

These  are  used  to  determine  the  equivalent  standard  deviation,  Sc,  in  column  10 
and  the  equivalent  skew  coefficient,  G^,  in  column  11.  The  standard  deviations 
and  skew  coefficients  for  the  hurricane  and  non-hurricane  frequency  curves  are: 
0.1330,  -0.8,  0.87  and  0,  respectively.  The  Wilson-Hil ferty  (1931)  approxima¬ 
tion  shown  in  Equation  4.3,  is  used  with  the  equivalent  skew  coefficient  of  the 
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Combined  curve  Gg  in  column  11,  and  the  normal  deviate,  Kp  in  column  7  to 
obtain  the  equivalent  Pearson  deviate  Kg  p  shown  in  column  12. 

KG,P  =  {[KP  ‘  V6)(GC/6)  +  1]3  '  1}  (2/6c)  (4.3) 

This  Pearson  deviate  is  then  used  to  develop  and  in  columns  16  and 
17  using  equations  9-4  through  9-6  in  Bulletin  17B  as  shown  below  for  the 
upper  and  lower  confidence  limits  respectively. 

zc 

Col .  13  =  3  =  1  -  yy 

2  4 

Col.  14-  b-  (KG>p)2  -  £ 

Col.  15  =  J  (Kgp)2  -  ab 


(4.4) 

(4.5) 

(4.6) 


(4.7) 


The  terms  Ep  c  and  Ep  c  in  columns  18  and  19  are  found  using  equations  3.8  and 
3.9.  These  are  then  used  with  the  standard  deviation  in  column  10  to  determine 
the  upper  and  lower  confidence  limits  Qp  c  and  Qp  c  with  equation  3.5  as 
shown  in  columns  20  and  21.  These  curves  are  drawn  as  a  dashed  line  in  Figure 
4.16,  revealing  the  irregular  shape  of  the  lower  confidence.  This  is  caused 
by  the  large  change  in  standard  deviation  near  the  intersection  of  the  hurricane 
and  non-hurricane  frequency  curves. 
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Method  2 


The  second  method  for  developing  confidence  limits  for  the  combined* 
population  frequency  curve  requires  that  separate  confidence  limits  be  cal* 
culated  for  the  hurricane  and  non*hurr1cane  frequency  curves.  Figure  4.17 
illustrates  these  computations.  The  5  and  95  percent  confidence  limits* 
shown  In  columns  3,  4,  6,  and  7*  are  calculated  using  Equations  9*4  through 
9*6  in  Bulletin  17B  (U.S.  Mater  Resources  Council,  1981)  for  the  discharges 
contained  In  column  1.  These  confidence  limits  are  drawn  on  Figure  4.18  and 
the  probabilities  associated  with  the  5  and  95  percent  confidence  limits 
(columns  9,  10,  12,  and  13)  are  selected  that  correspond  to  the  discharge 
shown  In  column  8.  Finally  Equation  3.4  Is  used  to  combine  the  hurricane 
and  non-hurricane  curves  yielding  the  combined  confidence  limit  curves  which 
are  tabulated  in  columns  12  and  14  and  shown  In  Figure  4.18. 

For  this  station  the  confidence  limits  calculated  using  Method  2  seem 
more  reasonable  than  the  one  calculated  using  Method  1. 
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Appendix  I 


DERIVATION  OF  EQUATIONS  3.3  &  3.4 


Fundamental  equation  representing  the  combined  probability  of  two  independent 
probabilities  is: 


x  P 


2 


This  represents,  however,  the  probability  that  both  P,  and  P2  will  occur. 
For  the  application  of  combining  two  frequency  curves,  the  question  is 
whether  one  or  the  other  event  will  occur. 


The  probability  of  occurrence  of  one  or  another  event  is  equivalent  to 
[1  -  (Probability  of  both  not  occurring-)] 

Probability  of  "non-occurrence"  for  P1  is  1  -  P-|  and  for  P2  is  1  -  P2 

Combined  probability  of  both  not  occurring  *  (1  -  P.j)(l  -  P2) 

Probability  of  either/or  occurring  =  1  -  (1  -  P-j )(1  -  P2) 

This  is  the  form  of  equation  3.3. 


Equation  3.4  is  found  by  algebra: 
Pc  =  1  -  (1  -  P^O  -  P2) 

.  1  -  (1  .  P1  .  P2  +  P,  p2) 


P2- 


P1  P2 
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